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Cowie (1977) and Cox and Anderson (1982, CA) have pointed out 
that when non-Coulomb processes heat the electrons In the shock 
front of a blast wave, strong thermal conduction at early times 
leads to a nearly Isothermal electron temperature distribution, 
but that at these same early times [t « t^^ = 5000 
yr(E^/10^^ ergs)^^^^(l cm"^/n^)^^^ for uni form external density] 

Coulomb collisions are slow to exchange energy between electrons 
and Ions and the Ion fluid remains essentially adiabatic. 

In the present paper, we apply methods developed In CA to 
this problem, utilizing an approximation scheme due to Kahn (1975). 

In so doing, we find the early self similtu' blast wave structure 
with adiabatic Ions, Isothermal electrons and equal electron and 
Ion temperatures at the shock. In subsequent papers, we shall pre- 
sent the time development of the system away from this early time 
form as Coulomb interactions become important and thennal conduction 
fades in significance (Cox and Edgar, 1982, Paper II). In addition, 
we shall complete the project begun In CA and Cox and Franco (1981, 

CF) presenting the X-ray appearance of a large set of possible 
explosions viewed from within, considering the possibility that the 
soft X-ray background derives from our being within such an'explo- 

_3 

slon. (The appropriate parameters, from CA, appear to be n^ 0.004 cm , 

5 50 

-V 100 pc, t 'V 10 years, for 5 x 10 ergs.) 

Apart from Its role in furthering the project mentioned, the 
present results will help to erase a certain parochialism sometimes 
present In Important discussions about the nature of shock fronts 
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in the presence of strong thermal conduction. In addition it pro- 
vides analytic results against which two-fluid numerical hydro- 
dynamic flow codes can be checked. 

Note that CA considered the effects of non-zero ambient 
pressure on a blast wave. Those effects are neglected in the 
present treatment, although for the soft X-ray background model, 
the evolving remnant is in an extremely difficult regime for which 
the external rressure, thermal conduction, and both Coulomb and 
non-Coulomb heating of electrons are all potentially of comparable 
importance. 

Turning to the problem at hand, the shock wave jump conditions 
for zero external pressure, generalized in CA to include thermal 
conduction flux arriving at the front are 



Pc = 
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where p, p, and u are pressure, density, and mass velocity, respectively, 
and the subscripts s and o denote post-shock and pre-shock values. 
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In addition is the shock velocity and x is the normalized density. 
Taking as the shock radius and assuming power law density be- 
havior, p « R,”“, the derivatives of the first three lines of 
O S 

equation (1) with respect to can be written 



Uj* = Xj*/(Xj-l) + Vj* (2) 

Ps* • P»* * 

where the asterisk is used to abbreviate the logarithmic derivative: 
fg* = dtnf^/dJinRj.. 

In the present context, with a self-similar evolution, x^ is 

3 

constant and x^* = 0. In addition, p^R^ is constant so p^* = -3. 

We shall retain the x^ dependence of the equations in anticipation 
of Paper II but set p^* = -3, as this will assumed to remain nearly 
valid even as the system evolves away from the early time asymptote 
explored here. The second two lines of (2) can then be rewritten 

X * 

u * = [u-3 + t^]/2 

[-3 - ^]/2 
s 

When the equations of motion are expanded about the shock 
position as in CA, the relations between the post-shock derivatives 
of the variables (still assuming p^* = -3) are found to be: 


(3) 
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where <>- - ~2 P) divergence of the thermal conduction 

flux. As can be seen by comparing equations (1) and (6), the jump 
conditions are directly related to the thermal conduction flux 
arriving at. the shock, but the post-shock derivatives depend on 
the divergence of F just behind the shock. Using the jump con- 
ditions, the left hand side of (6) can be rewritten 
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Knowledge of the parameter ts essential to knowlecige of 


t’le post-shock derivatives, but cannot itself be known without 


evaluating /R tdR over the entire structure. 
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The conservation of energy equation 
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can be separated Into Its two fluid components 
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where the energy exchange rate between Ions and electrons Is 
(Spitzer, 1962-, itoh, 1978; CA) 


h > akn^(T. - = ftf 
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Here k is Boltzmann's constant, n " a tnA/153 cgs. 


inA = An[1.2 X 10^(T^^^Tg/n^^^)],nj^g 


0.1 n„ was assumed, T. is 

n I 


the ion temperature, the electron temperature, and T = 

(1.1 +1.2 Tg)/2.3 Is the average temperature. Similarly, = 

nkTi* Pg * (1.2/l.l)nkTg, p = (2.3/1 .1 )nkT. 

These separated equations can be manipulated Into a variety 

1 

of useful forms. For example, defining g » T^/T, 
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The second term on the right side is the Coulomb heating and is negligible 
for times small compared to the equilibration time mentioned earlier 
(see CF for values when w 0), unless the first term is also small. 

The first term is conduction cooling and will be small at early times 
either if the electron temperature is not equilibrated behind the shock 
or thermal c;>nduction is thoroughly quenched by a magnetic field. 

Itoh's (1978) useful approximation to the electron temperature distributions 
for adiabatic models follows from setting & » 0 and recognizing that 
n/T^/^ constant for a parcel. The g's from the RHS of (11) 

can be included in a total time derivative on the IMS which is then 
constant. In the present study, a more useful result is 
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together with a slightly altered version of (8); 
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In this form, two separate variables are identified, one of which is 
sensitive only to Coulomb heating, and one of which is sensitive only 

to i. 


Since we are assuming that at the shock front, = T, due 

to non-Coulomb heating processes, and are interested in the similarity 
solution for small times, we will neglect the RHS of equation (12). 

Near the center of the structure, g will be very small so that it would 
appear that this neglect is not entirely justified. But in this case 
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equation (12) implies that 0/0t[tn(T^^^/n)] « n, which vanishes for 

small R, so the appv'oximation holds. 

For t « tg^, then, a particular parcel of gas has two constants 

of the motion which allow it to be located within the structure and partly 
constrain its properties. One of these constants is M(R), the mass 
enclosed by the radius through the parcel. The other is the dif- 
ferentiated quantity in (12). 

A convenient normalized mass coordinate is p * M(R)/M(R^). When 
a parcel of current mass coordinate v was first shocked, its temperature 

was T^(R_)/y where T^(R^) is the current value of the post-shock 

temperature (at p ■ 1). Since M(R^) « R^'’"“, the parcel was originally 

1 

located and shocked at R^(;i) p'’”“ R^. Its density at that time was 

X (R.) n (R.) * X (R,)(R^/R 4 )“ n (R,). Since we are considering only 
the early self similar form, x^(R^) = ' ’^s* ' 

x^n^(R^)p’“^^^”“^. In addition, the value of g just inside the shock 

was 1, so that the second constant of the motion can be written in terms 
of p and current values of T. and n , 
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We are furthermore exploring an example for which thermal conduction 
is sufficiently rapid that everywhere, making g » Sub- 
stituting this into (14), recalling x ■ solving for T, 


T 


Li 4. Li 

2.3 2.3 
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where x is the current normalized density at the location of the mass 
element with present mass coordinate p. 

By conservation of mass, au/3r ■ (3-ai)r x where r ■ R/R^. M^th 
this result, the local value of « ran br derived fron >quation (13). 

A number of steps are involved in the algebra, including expansion of 

D 3 d 2 

^ ^ ^ result u « [r - y/(r x)] v^ which follows 

from mass conservation. The net result after several fortuitous can- 
cellations is (with x' = 8x/3r) 


^ ^ \ r^x 




(16) 


1 2 

A self-consistent solution requires ^ r [-y]dr * 1 since is the 

integral of the divergence i. The first term in (16) is easily integrated 

1 2 -1 * 
since ^ r xdr - (3-w) but the second term still requires a complete 

knowledge of the structure. 

Evaluating equation (16) at r » u = 1 , x'/x^ * o, however, provides 
another relation between the post shock derivatives 
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wh1ch» when combined with (5), (6), and (7) with ■ 0 can be solved 
for a, B. 0 and R.t /F as functions of x alone. (There are four 
parameters and fou<^ equations, one each from conservation of mass and 
momentum, and two from the separated energy equations.) There Is then 
a one parameter family of edge conditions only one of which will 
correspond to a complete solution satisfying the Integrated flux condition. 

Thus far no approximations regarding the structure have been made. 
Three possible approaches can be taken at this point One Is to put tr.e 
exact equations of motion and boundary conditions on a machine and ask 
It for the answer. A second is to df'velop a method For findinq a 
reasonable approximation to ^•'“om which all of the post shock 

parameters and their first derivatives would follow automatically. A 
third Is to apply Kahn's (1975) technique to approximate the entire 
structure, given x^ and the parameters it determines, and then to 
perform the flux Integral as a function of x^, choosing the value of x^ 
which effects the flux normalization. 

We have carried out the last two of these approaches, the second 
because the estimate helps to understand the results, the third because 
(1) Kahn's (1975) approximation technique has been shown to produce 
extremely accurate results in the adiabatic problem, (see appendix to CF); 
(ii) because the structure can be described analytically; and (1i1) 
because this scheme Is used in Paper II to find the later behavior so 
what we actually require for Paper II Is the asymptotic behavior of the 
approximation. The truth, of course. Is that we have done the problem 
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we felt competent to solve, confident that the answers would be close 
to correct and could be presented conveniently. 

In attempting a first approximation to we note that 

all electrons have the same temperature at any given time, Implying 
that they all lose energy at an equal rate. If all losses were via 
thermal conduction, i/p would be a constant, Flux normal Iza 

tion Integral can be written 


or 
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SO we expect the quantity in square brackets to be of order unity, and 

we have found a reasonable est1me<.e for Electrons also lose 

energy by expansion, however, and these losses are most severe near the 

edge. Thus the edge electrons necessarily have i/p less than the mass 

average In order to have the same total loss rate. 

In order to obtain an estimate of the size of this effect, an analysis 

similar in spirit to the ore presented here was carried out for a totally 

isotheniial blast wave, the solution to which has been given by Sollinger, 

‘s t 

Rappaporc, and Buff (1975, SRB). Assuming that — ■ led to x^“ 2.457 

whereas SRB quote x ■ 2.378, a difference of only 3.3X. Using x ■ 2.378 

to solve for the ratio of edge to mass-averaged l/p, R i /[F x (3-w)l, 

ss ss 
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the result was 0.862. Thus the greater Importance of expansion losses 
near the edge leads to an only 14t lower than average. 

Returning to the present Investigations, when all the substitu- 
tions are made to utilize the mass averaged t/p, the results can be 
expressed as 




X 


36[(6x^ -llx,-4) - «x^(4x^-7)] 

(3-<o)(4-x:)x i9rx-:T^Tn — 


(19) 
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where some of the peculiar mmibers result from the inclusion of lOX 
helium by number. For u = 0, the ratio equals 1 for = 3.3456 or 1.5176. 
As u (explosions In cavities have negative u), these values 

approach = 2.3972 and 1.779 respectively, fur the totally 
isothermal case, the equation corresponding to (19) differs only In 
that (91x^-160) In the denominator is replaced by 36{x^-2) in which 
case the two solutions for u » 0 and ratio 1 are x^ ■ 2.457 and 1. 

We have not explored the significance, if any, of the lower compression 
possibilities, for example, whether or not they could correspond to 
thermal waves. It is clear, however, that the higher compression 

I 

solutions of each pair are the analogs of the x ■ 2.378 solution of 


SRB. 


Notice that there Is no single value of x^, u,., or other post 
shock parameter which Is always appropriate to shocks with Infinite 
thermal conductivity behind them. The sh.-.k compression depends on 
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the complete structure, here summarized by the u dependence. Notice 
too that as pointed out by Cowle (1977) the SRB result exaggerates 
the deviations expected from the usual adiabatic structure If the 
Ions remain adiabatic. 

As previously intimated, the results can be Improved upon 
considerably and the complete structure given by resorting to Kahn's 
(1975) approximation scheme for the structure, generalized as In 
CA and CF. We first solve for o(x^) using equations (17). (5), (6), 
and (7) with x^* * 0 (again with ■ 0.1); 

wx 

,(x -1)[(20.?5)x - 13.8] - -a [(17.9)x -31.7], 
a - 9.V8 ^ 

A second parameter required Is 


X (3+a) - (3-«)x ^ 

x^-5/2 

which was derived by matching the form of the structure to the post 
shock density slope, a. The resulting structure Is given by 


,(r) . tr^'^xp ((Xj-|) 1^1)]’-“ 
x(r) . [| * (Xj-f)r'’] 4 «(>■) 


( 22 ) 


with T(r) provided by (IS). This form automatically conserves mass, 
provides the correct post shock values of all parameters (when 
supplemented by the jimip conditions) and gives the correct post shock 
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derivatives. It also correctly matches logarithmic derivatives in 
the limit of small r. All that is required is to find the correct 
value of Xj. 

In this approximation, the normalization integral for the flux 
becomes equivalent to Q(x^, u) = 1 where 
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1.1 1 

2.3 4-Xj 


, 3(1-0.) 
' 3 - 0 . 
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These integrals have been evaluated numerically for o. = 0, -2, and -4, 
varying x^ until the Q = 1 condition was satisfied. The resulting 
values of x^, o, q, and given in Table I. 

Several comparisons between the two levels of approximation are 
in order before considering the complete structures associated with 
these results. First we note that values of x^ obtained by setting 
equation (19) equal to 1 are 3.3, 3.1, and 2.6 per cent larger (for 
0 . = 0, -2, and -4 respectively) than those in Table I, the same degree 
of accuracy obtained for the isothermal case. Secondly, the values of 
in Table I are very well approximated by Xj(2.53-to) while the 
mass average is x^(3-u). Therefore, cases which we have not presented 
should be closely approximated by setting equation (19) equal to 


(2 . 53-(i.)/(3-u) . 
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Finally, the density, pressure, temperature, and Rgl/[F^x] 
(normalized t/p) structures are shown graphically for u » 0, -2, 
and -4 In Figures 1 through 3. They bear a strong resemblance to 
the adiabatic structures. The central pressure plateaus are slightly 
lower (p/p^ » 0.26, 0.35, 0.39 for u * 0, -2, -4 respectively) than 
when the same approximation Is applied to the adiabatic case 
(p/Pj “ 0.31. 0.38, 0.40). The edge densities are lower 
(x^ • 3.24, 2.91, 2.77) than the x^ ■ 4 for the adiabatic case, but 
the density gradients are less so that the Interior densities are 
higher hy factors of 1 34, 1.41 , 1,44 for 0, , 4. for w - 0, 

the crossover occurs near R/R^ * 0.9, x » 1 .5. Dividing the pressures 
by the densities Implies that the edge temperatures are higher by 
factors of 1.23, 1.37, and 1.44 for the same p^, and p^, but the 
central temperatures (local mean, T) are lower by factors of 0.62, 

0.66, and 0.68 relative to the adiabatic structures. For given 
total energy, R , and p , we shall see that p differs by at most a 
few per cent from the adiabatic case so these ratios are altered 
very little. 

In order to know the explosion energy represented by a, particular 
solution. It is necessary to Integrate the thermal and kinetic energies 
over the structure and sum them. The thermal energy Is 
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where Ij ■ 3 ^ r^(p/p^)dr, and 
Similarly, the ki'^etlc energy is 


4irR^^/3. 


E„ = 4irR^p(u^/2)dR 
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where = (3-w) ^ r^x(u/u^)^dr. 

The results of numerical integrations of Ij and I|^ for o» ■ 0, -2, 
and -4 are provided in Table II. Table II also lists E^/(3V^p^/2) for 
these same value r w (where E = Et + fi,) . as well as the corresponding 
values from CF for the completely adiabatic structures. • 

From these results and the jump condition for p^, the time 
development can be found. For example with w = 0, = constant, x^ = 3.2383, 

(x -D/x. » 0.6912 p V ^ 3V D /2 - 4.343 R. Vv ^ 

= E^/0.6430 (from Table II), or - 2.793 » 0.447 

(since for u * 0, = 0.4 R^/t). Thus, R^ » (2.238 E^t^/p^)'^^. 

The adiabatic result is similar, but with 2.238 replaced by 2.03. 

Thus strong electron thermal conduction in the *• 0 structure increases 


Rs and Vj by only about 2% over an adiabatic system. 


r\h 

In the approach taken in Paper II, It Is necessary to make an 

a priori choice for the behavior of p^*, after which a self-consistent 

solution Is found for the equations for motion (again assuming the 

Kahnian density distribution). Any fault of the choice for p^* then 

appears as a time variation of the systemic total energy. 

We note that as the system evolves away from the early time 

asymptote, It Is evolving toward equilibrated electron and ion 

temperatures and toward negligibility of the thermal conductivity 

(e.g. CA). Hence (until radiative cooling becomes Important), the 

late time asymptote Is Just the adiabatic structure discussed In CF. 

And the fact that the total energy integrals of Table II show very 

similar values of E /(3V,p,/2) for the two asymptotes guarantees that 

0 s s 

on average, p^* » ntnp^/atnR^ a -3 during the transition. In addition, 
major fluctuations in p^* seem unlikely since so little readjustment 
Is required in the density distribution. (Fluctuations in p^* would 
provide transient alterations In the pressure gradient, which In turn 
cause transient accelerations which rearrange the mass distribution.) 

The only possible difficulty may arise from the suddenness with which 
the transition Is made. Thus taking p^* « -3 throughout the transition 
Is a reasonable first step, although the transition structure must then 
be examined for nonconservation of energy. 

We summarize our results by way of example. For a particular 

appl Icatlor, t , u, and p for a particular value of R will be known. 

o o 

If 0 ) Is 0, -2. or -4, E /[3V o /2] • E /[2i.R,^/p J Is found from Table II, 

O S S 0 S S 

which then Is solved for 1^®"" other negative values of u, the 
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results can be Interpolated to reasonable accuracy. Similarly, values 
of x^. a, and q can be taken from Table I and substituted Into 
equations (22) and (15) to find the normalized density, temperature, 
pressure [from p ® (2.3/1 .1 )nkT], and mass coordinate distributions. 

The jump conditions, (1), together with Pj(R 5 )» provide 

v^(Rs) which can be Integrated easily to find Rg(t). Finally, the 
electron temperature Is everywhere 

Again, for values of u other than those In Table I, can be 
found by setting equation (19) equal to (2.53-aj)/(3-a)) , after which x 
and q follow from (20) and (?1). Tboreaftor solution proceods as before. 


The authors wish to thank Josd Franco for carefully reviewing 
much of the algebra underlying this work. This research was supported 
in part by NASA Grant NGL 50-002-044 at the University of Wisconsin. 
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TABLE I 


U) 

X 

s 

a 

q 

Rs4/Fs 

0 

3.2383 

8.6044 

8.2879 

8.176 

-2 

2.9121 

13.8907 

16.4663 

13.212 

-4 

2.7721 

18.9004 

25.4258 

18.097 
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TABLE II 


Present Cese Adiabdtic (CF) 



•t 


Energy® 

•t 

'k 

Energy* 

0 

0.474 

0.735 

0.643 

0.470 

0.740 

0.655 

-2 

0.473 

0.800 

0.578 

0.465 

0.773 

0.581 

-4 

0.434 

0.844 

0.511 

0.464 

0.790 

0.549 


Eo/[3V^p^/2] = 1-f + (x^-l)I|^/[x^(3-u))] 
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FIGURE CAPTIONS 


Fig. 1. — Self-similar structi're of a Llastwave with p ■ constant, 

O 

adiabatic ions, isothermal electrons, and at the shock. 

Xj = 3.2383, 

a) p/p. and x/x versus R/R 

b) 




p/p^ 



as functions of the relative mass coordinate u. 


Fig. 2 . -- Same as Fig. 1 but with p « R, 

0 s 

.4 

Fig. 3. — Same as Fig. 1 but with • 


2 . 9121 . 

2.7721. 
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